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Abstract 

We consider some basic principles of fluid-induced lubrication at soft interfaces. In particular, 
we show how the presence of a soft substrate leads to an increase in the physical separation 
between surfaces sliding past each other. By considering the model problem of a symmetric non- 
conforming contact moving tangentially to a thin elastic layer, we determine the normal force in 
the small and large deflection limit, and show that there is an optimal combination of material 
and geometric properties which maximizes the normal force. Our results can be generalized to 
a variety of other geometries which show the same qualitative behavior. Thus they are relevant 
in the elastohydro dynamic lubrication of soft elastic and poroelastic gels and shells, and in the 
context of bio-lubrication in cartilaginous joints. 
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Lubrication between two contacting surfaces serves to prevent adhesion and wear, and 
to reduce friction. The presence of an intercalating "lubricating" fluid aids both, but gives 
rise to large hydrodynamic pressures in the narrow gap separating the surfaces and can 
thus lead to deformations of the surfaces themselves. For stiff materials such as metals 
the pressure required for noticeable deformations is very large (~ lGPa) and under these 
conditions the lubricating fluid will exhibit non- Newtonian properties . However, if these 
surfaces are soft, as in the case of gels and thin shells, elastohydrodynamic effects can become 
important when the fluid is still Newtonian since the pressure required to displace the surface 
is appreciably less. This type of situation is also common in mammalian joints where the 
synovial fluid serves as the lubricant between the soft thin cartilaginous layers which coat 
the much stiffer bones. Motivated by these observations, in this letter we consider the 
coupling between fluid flow and elastic deformation in confined geometries that are common 
in lubrication problems. 

As a prelude to our discussion we consider the steady motion of a cylinder of radius R 
completely immersed in fluid and moving with a velocity V, with its center at height ho + R 
above a rigid surface (see Fig. 1). The dynamics of the fluid of viscosity /i, and density p 
are described using the Navier-Stokes equations 

p(d t v + v ■ Vv) = /iV 2 v - Vp, (1) 

V-v = 0, (2) 

where v is the 2-D velocity field (u, w) and p is the pressure. Comparing the ratio between 



the inertial and viscous forces in the narrow gap having a contact length / ~ \/R1iq [14] , we 
find the gap Reynolds number Reg = ^jjp ~ p ^i° /2 <C = Re, the nominal Reynolds 
number. The gap Reynolds number is small since ho <C R and we can safely neglect the 
inertial terms and use the Stokes' equations (and the lubrication approximation thereof 
P]) to describe the hydrodynamics. The temporal reversibility associated with the Stokes 
equations and the symmetry of the parabolic contact leads to the conclusion that there can 
be no normal force due to the horizontal motion of the cylinder. However, if there is a thin 
soft elastic layer on either the cylinder or the wall, the deformation of the layer breaks the 
contact symmetry and leads to a normal force. This then leads to an enhanced physical 
separation and a reduced shear so that it may be a likely cause for the low wear properties 
of cartilaginous joints. 




FIG. 1: A rigid cylinder moves at a velocity V a distance Hq above a rigid substrate coated with 
an elastic layer of thickness Hi. Hi, ho <C s/HqR = I. We illustrate the steps of the perturbation 
analysis: (b) an antisymmetric pressure distribution pushes down on the gel in front of and pulls 
the gel up behind the cylinder; (c) the fore-aft gap profile symmetry is broken; (d) the new pressure 
field produces a normal force, (a) and (b) correspond to an undeformed substrate, while (c) and 
(d) correspond to solutions of Q), © and (JT5J) for rj = ^ = 10. 

Continuing our analysis in the context of a cylinder moving along a planar wall, we take 
the direction to be parallel to the wall in the direction of motion of the cylinder and the 
z— direction to be perpendicular to the wall; p is the fluid pressure; h is the distance between 
the solid surfaces. Guided by lubrication theory 3] we use the following scalings 
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to reduce (JTJ and (0) to 



d x P = d zz U, d z P = 0, 
dxU + dzW = 0. 



(3) 

(4) 
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We consider steady motion in the reference frame of the cylinder, so that the boundary 



conditions are 



U(X,0) = -1, U(X,H) = 0, 

W(X,Q) = W(X,H) = Q, (6) 

Integrating (I4l5l(il) leads to the dimensionless Reynolds equation Q]: 

= d x (6H + H 3 d x P). (7) 

Since the gap pressure is much larger than the ambient pressure, we may approximate the 
boundary conditions on the pressure field as 

P(oo) = P(-oo) = 0. (8) 

Next, we consider the deformation of the elastic layer of thickness Hi that rests on a rigid 
support. Balance of stresses in the solid leads to 

V • a = 0, (9) 

with the stress given by 

a = G(Vu + Vu T ) + AV • u I, (10) 

where u = (u x ,u z ) is the displacement field and G and A are the Lame constants for the 
solid, which is assumed to be isotropic and linearly elastic. To calculate the increase in gap 
thickness H(x) we use the analog of the lubrication approximation in the solid layer 4j. The 



length scale in the z— direction is Hi and the length scale in the x— direction is \/hoR. We 
take the thickness of the solid layer to be small compared to the thickness of the contact 



zone, yh^R 3> Hi, and consider a compressible elastic material, G ~ A, to find the vertical 
force balance: d zz u z = 0. The boundary condition at the solid-fluid interface is er • n = — pn, 
so that (2G + X)d z u z (x, 0) = —p(x). Using the zero displacement condition at the interface 
between the soft and rigid solid, u z (x, —Hi) = leads to the following expression for the 
displacement of the surface 

Hip(x) 

Then the dimensionless version of the gap thickness, h = h + |^ — u z (x, 0), is 

H(X) = 1 + X 2 + V P(X), (12) 
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where ri = Ah/ Tin = -^0EMm}L j s the dimensionless parameter governing the size of the 
deflection. Inspired by the some recent experiments p| in a similar geometry, we consider 
a cylinder of radius R = 10 cm coated with a rubber layer (Hi = 0.1 cm, G = 1 MPa) 
moving through water (p = 1 mPa-s, V = 1 cm/s, h = 10~ 3 cm). Then rj = 10 -2 << 1, so 
that we may use the perturbation expansion P = P + f]P\, where Pn is the anti-symmetric 
pressure distribution corresponding to an undeformed layer, and Pi is the symmetric pressure 
perturbation induced by elastic deformation. Substituting (fT2|) into (J7|) leads to the following 
equations for Pn, Pi: 

if: d x [Q(l + X 2 ) + (l + X 2 ) 3 d x P } =0, (13) 
ri 1 : d x [QP + 3(1 + X 2 ) 2 P d x P + 

(l + X 2 fd x P 1 } = 0, (14) 

subject to the boundary conditions Pn(oo) = P (— oo) = Pi(oo) = Pi(— oo) = 0. Solving 
(JED) and (JUJ) yields 

^" (l+X2)2 +r/ 5 (l+X2)5- ^ 

Then the normal force is 

3tt 



PdX = — V} (16) 

~° n 

In dimensional terms, F = 3 ^ 27T M 7 V 9 HlR 1 , whose scaling matches the result reported in ml, 
but with a different pre-factor. When rj is not small, we solve jZj),© and (fT2"|) numerically. 
Figure El shows that as rj increases the mean gap increases and its profile becomes asymmetric, 
resembling the profile of a rigid slider bearing, a configuration well known to generate lift 
forces 3]. In addition, this increase in the gap size causes the peak pressure to decrease 

since p ~ ^ V f/i ■ These two competing effects produce a maximum lift force when t] = 2.06. 

V 

The physical basis for the previous arguments can be more easily understood using scaling 
and therefore allows us to generalize these results to a variety of configurations involving 
lubrication of soft contacts (Fig. 3; Table 1). Balancing the pressure gradient in the gap 
with the viscous stresses yields 

p pV pVR 1 ' 2 

l ~ h 2 ?~ h 3/2 > V) 

Substituting h = h + Ah, with Ah <C h , we find that the lubrication pressure is 

fiVVRh~o pVR 1/2 n , Ah Ah 

P~ 77 , AM2 777- (1 + 7— ) =P0 + 7— Pi- (18) 

(h + Ah) 2 ^/2 ho h(j 
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FIG. 2: (a) The dimensionless pressure distribution, P, for several values of rj, a measure of 
the deflection of the elastic layer compared to the initial separation, (b) The dimensionless gap 
thickness profile, H = 1 + X 2 + rjP. The gap thickness and asymmetry increase with rj, while 
the maximum value of the pressure decreases, (c) For small r] asymptotic analysis predicts a 
dimensionless lift force F = ^77, which matches the numerical solution, (d) F has a maximum at 
rj = 2.06 as a result of the competition between symmetry breaking and decreasing pressure. 



Here p does not contribute to the lift for the reasons outlined earlier, so that the the lift 
on the cylinder per unit length is 



Ah fxVR 

'Pi 1 L2~ Ah i 

fin 



h 



(19) 



where Ah is determined by the solution of the elasticity problem. 
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FIG. 3: Schematic diagrams of two configurations considered on the level of scaling: (a) a soft solid 
coats a rigid solid where Hi 3> \/hoR, i.e. the layer thickness is larger than the length scale of the 
hydrodynamic interaction; (b) the cylinder is replaced by a cylindrical shell. 



TABLE I: Summary of scaling results for small surface deflections. 



Geometry 


Material 


Surface displacement 


Lift force/unit length 


Thin layer 


Compressible elastic solid 


fiV HiR 1 / 2 

G ,3/2 


^V 2 H1R 3 / 2 
G h 7/2 


Thin layer 


Incompressible elastic solid 




^V 2 HfR 1 / 2 


G h 5 n /2 RV 2 


G h 9/2 


Thin layer 


Poroelastic solid 


\iV HiR 1 / 2 

Geff h */* 


M 2 V 2 H t R 3 / 2 

Geff h y 2 


Thick layer 


Elastic solid 


iiV R 
G h 


V 2 V 2 R 2 

g 4 


Cylindrical Shell 


Elastic solid 


txV R 7 / 2 


ll 2 V 2 R?/ 2 




g Kh y 2 



For a thin compressible layer, the case treated above, the normal strain is ~ ^ ~ 
~^TT- Therefore, 

liVHiR 1 ' 2 „ n 2 V 2 HiR 3/2 , s 

A/i ~ F ~ (20) 

G h^ 2 ' G h 7/2 ' 

In sharp contrast, a thin incompressible layer will deform via shear with an effec- 
tive shear strain ^ ~ ^[15]. Balancing the elastic energy J G( R 1 ^ Ah ) 2 dV ~ 
G( - ' 1 \ Afc ) 2 i7;v / Rho with the work done by the pressure p Ah^Rh yields Ah in terms 
of p Q . Then, JEJ and (dHJ) give 

A thick layer (H\ ^> y/h R) may be treated as an elastic half space. The strain scales 
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as Jf^ R and remains appreciable in a region of size h R. Balancing the elastic energy 
/ G{^=) 2 dV ~ G(-^=) 2 Rh with the work done by the pressure p Ahy/Rh~ yields Ah 
in terms of po- Then, (fTTjl and (fTT?|) give 

. pV R fi 2 V 2 R 2 

A ^^V F ~—J4 (22) 

Finally, we consider the case of a cylindrical shell, of radius R and thickness h s , moving 
over a rigid substrate. Since the shell is thin it can be easily deformed via cylindrical bending 
without stretching. The bending strain is hs ^ h so that the elastic energy J G( hs ^ h ) 2 dA ~ 
Ghs R? h ■ Balancing this with the work done by the pressure Po(\/Rho) 2 ^ yields Ah in terms 
of p . Then, flT7|) and (JUJ) give 

F ~—l^ (23) 

We note that the above scalings for cylindrical contacts can be trivially generalized to 
spherical contacts for the case of small deformations, but space precludes us from discussing 
these in detail. 

We conclude with a discussion of how our results may be applied to the lubrication of 
cartilaginous joints lllsl, where a thin layer of a fluid-filled gel, the cartilage, coats the stiff 
bones and mediates the contact between them. Here, electrostatic effects prevent physical 
contact of the surfaces under high static normal loads, while elastohydrodynamic effects 
could enhance separation and thus reduce wear. Inspired by the treatment of cartilage using 
poroelasticity , the continuum description of a material composed of an elastic solid 
skeleton and an interstitial fluid , we treat the cartilage layer as an isotropic poroelastic 
material |l6[. The gel can then be described by its fluid volume fraction a ~ 0(1), drained 
shear modulus G and drained bulk modulus K ~ G, thickness Hi, permeability k, and 
interstitial fluid viscosity fi. Using dimensional reasoning, we can construct a poroelastic 
time scale 



(24) 



which characterizes the time for the diffusion of stress over the layer thickness Hi due to 
fluid flow. Then the response of the gel is governed by the relative size of r p to the time scale 
of the motion, r ~ l/V = y/h R/V. If r ^> r p , the motion is so slow that the interstitial 
fluid plays no role in supporting the load. If r ~ r p , the fluid supports some of the load 
transiently, thereby stiffening the gel. Finally, if r t p the response of the gel will depend 
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on the size of the Stokes' length l s ~ \J rpj p. If l s ~ Hi, there is no relative motion between 
the fluid and the solid and the gel behaves as an incompressible elastic solid ljj, ll2j, with 
shear modulus G. From (|2U|) and (J2TJ) we see that the effective modulus is 

G eff ~ po^/M ~ -^G ~ -j^fG. (25) 

To find the scale of the deflection and lift force we use the same scaling analysis as for a 
thin compressible elastic layer but replace G with G e ff(r), so that (J2DJ yields 

Ah ~ — F ~ — 5 (26) 

where G e ff E [G, j^G}. Inserting characteristic values V = 1 cm/sec, G = 10 7 g/sec 2 cm, 
Hi = 0.1 cm, R = 1 cm, /i = 10~ 4 cm, - = 10 -13 cm 3 sec/g shows that r <^ r p , but since 
I ~ ff; significant material stiffening is prevented. Consequently, the effective modulus is G 
and the scale of the deflection is 

yV HiR 1 / 2 , 
" ~ ~ (27) 

which suggests that joints could easily operate in a parameter regime that optimizes repulsive 
elastohydrodynamic effects. Although our estimates are based on non- conforming contact 
geometries; in real joints where conforming contacts are the norm, we expect a similar if not 
enhanced effect. 
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